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Abstract 



The paper presents equations determining the particle spin evolution in 
the post-Newtonian approximation in the problem of motion of two mass 
and spin possessing particles. The equations are derived with the Einstein- 
Infeld-Hoffmann method from the condition of metric tensor symmetry. The 
consideration uses the condition of coordinate harmonicity and the metric 
coincidence nearby the particles with expansions of the Kerr solution written 
in the harmonic coordinates. For gyroscopes on satellites Gravity Probe 
B, for example, the equations yield a deviation of the axis of revolution, 
which is, first, two times as small as that obtained by J.L.Anderson in paper 
gr-qc/0511093| and, second, the deviation is of opposite sense. From the 



equations it follows that the total angular momentum of the spin particle 
system, generally speaking, is not conserved, beginning even with the post- 
Newtonian approximation. The results are briefly discussed. 
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1. Introduction 



A system of two particles possessing masses and intrinsic angular momentums 
is considered. The description of motion of these particles within the general 
relativity requires two types of dynamic equations: the ones determining the 
motion of center of masses of particles and spin evolution for each of them. 

The equations of the first type can be derived either with Fock method 
or Einstein-Infeld-Hoffmann (EIH) method. As a rule, no disagreements in 
regard to these equations arise. 

The situation regarding the spin evolution describing equations is ab- 
solutely different. There are several versions of these equations for self- 
consistent motion of two spin particles (e.g., [l]-[4]) and plenty of versions 
of equations describing the spin dynamics of the trial spin particle (e.g., [5]- 
[17]). All the versions differ from each other, and it is far from being evident 
what of these versions is proper. 

We have solved the problem of self-consistent motion of two spin particles 
in the post-Newtonian approximation (PN^ ) with the EIH method. The 
distinctive features of our consideration are the following. 

1. All the operations provided for by the EIH method are performed 
using de Bonder condition for metric 



which coincides with the condition of coordinate harmonicity. 

2. In the approximation expressions of the EIH method the unknown 
coefficients involve expansions of the Kerr solution written in harmonic coor- 
dinates in [18], [19]. The Kerr solution expansion coefficients are taken from 



3. Dynamic equations for the spins result from the conditions of the 
metric symmetry in the approximations under consideration. 

The dynamic equations for particle spins that have been derived by us 
and the ones for trial spin particle which follow from them differ from all oth- 
ers in that they are consistent with the starting general relativity postulate 
of metric tensor symmetry. Ensuring the metric symmetry, we encounter an 
alternative: either the metric symmetry or the law of conservation of angular 

^By the PN approximation is meant an approximation, where corrections to the Newto- 
nian approximation begin to appear irrespective of the order of smallness of the corrections. 




(1) 



ref. [20]. 
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momentum. As a result of the choice made we arrive at nonconservation of 
the vector of total angular momentum of the particle system even in the PN 
approximation. This issue is interesting by itself, and we discuss it briefly. 
Besides, estimations for gyroscopes, the experiments with which were con- 
ducted on satelhte Gravity Probe B (GP-B), are presented and discussed. 

We guess that the reader is familiar with the EIH method, so we restrict 
ourselves to minimum elucidations in this regard. Also, cumbersome treat- 
ments on EIH technique are not included in this paper. All these treatments 
and calculation details will be presented in a separate paper. 

2. Notation 

This paper solves the equations of general relativity with zero energy-mo- 
mentum tensor. Like in refs. [21], [22], write the equations using rjafs, the 
metric tensor of Minkowski space. 

Raf3 - iv'^'R^.u) = 0. (2) 

Here Raf3 is Ricci tensor specified in the standard manner. In the back- 
ground space the Cartesian coordinates will be used, therefore 

Vafs = v"^ = dtag [-1,1, 1,1]. (3) 
Two types of quantities are used: hap, 'jais- These quantities are given 

by: 

The raising and lowering of index in hap and 7^/3 is made using tensors 
©. Thus, by h"^ is meant /i"^ = ri''^'ri^''h^y. 

A system of two particles is considered. In the background space we 
choose Cartesian coordinates [ct, Xk) in an arbitrary way. Denote the particle 
radius- vector coordinates by ik',Vk, respectively. Denote the length of the 
vector with coordinates Rk = i]k—^k by R. Denote the derivative with respect 
to ct either by point above the quantity (for example, ^k) or as subscript "0" 
separated with comma (for example, ^^^o = ^fc)- 

The masses of the first and second particles are denoted by M, fh, respec- 
tively. Quantities 
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, ^ GM Grh 
M = — —; m = 

c c 



are equal to halves of Schwarzschild radii of the particles. Here G is the 
gravitational constant, c is light speed. In what follows quantities M, m will 
be referred to as masses, although they have the length dimension. Write 
the axial vectors of the intrinsic angular momentums of the particles in form 



The Sk] Sk have the length dimension; we will call them the reduced vec- 
tors of intrinsic angular momentums or simply the spins. 

Equations © in the A'^-order are denoted as [00; A^], [0/c; A^], [mn; A^], etc. 
Hereinafter the coordinate conditions are denoted as [c.c; 0; A^], [c.c; fc; A^], 
etc. 

3. Smallness parameters 

In the EIH method the principal smallness parameter is 



where v is the characteristic relative particle velocity. Alongside the param- 
eter A, the system under discussion has four more dimensionless and a priori 
independent parameters: 



With large enough values of the radial variable r these parameters become 
small. However, the relative relationships among them are arbitrary in prin- 
ciple as they are independent of r and depend on the particle characteristics. 
In this paper the consideration will be under the following assumptions: 



McSk', mcsk- 



A = v/c 



(5) 



u = M/r; V = S/r; m/r; s/r. 



(6) 





(7) 



S 



s 



A. 



(8) 



r 



r 
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We will not dwell on elucidation of the physical meaning of these assump- 
tions; these elucidations are presented in many papers and are of standard 
nature. Note only that in many real situations assumptions (|7j), (jH)) are satis- 
fied not in full measure. Nevertheless, the dynamic equations derived under 
assumptions ((Tj), ^ can be used as a source of information about the motion 
of spin particles in these cases as well. 



4. Arrangement of orders of smallness 



Certain assumptions as to the arrangement of the orders of smallness of some 
terms appearing in functions 700, Jok, Imn should be made. The lowest-order 
approximation term construction is determined uniquely by the assortment 
of parameters The orders of smallness of the other terms are taken so 
that the EIH procedure be self-consistent. Our consideration will use the 
following arrangements of the orders of smallness: 

700 = 7 00 + 7 00 + 7 00 + 

2 4 5 

7ofc = 7ofc + 7ofc + 7ofc + \ (9) 

3 5 6 I ^ ' 

Imn ' mn + 7 mn + 7 mn + • • • ■ 

4 5 6 ) 

Typically, in the EIH method, the expansions of functions 700, 7ofe, 7mn are 
composed of terms of the same parity in A. In our case the results of [20] are 
taken into consideration, from which it follows that the expansions of func- 
tions 700,70^5 7mn should contain terms of all orders of smallness, beginning 

with 7 00, 7 ofc, 7 m„, respectively. 

454 

Recall some "techniques" used in the EIH method. 

• In the EIH method, in the differentiation of any function in x° its order 
of smallness increases by one as against the one that the function had before 
the differentiation. In the differentiation with respect to spatial coordinates 

the function order of smallness is preserved. 

• All the expressions must be symmetric about the replacement of the first 
particle by the second and the second by the first, that is about simultaneous 
replacements: 



M —>■ m; m M; S ^ s; s ^ S; 

^k] Vk ^k] Rk ~Rk- 



(10) 
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• With vanishing parameters m, s the expressions derived with the EIH 
method should be the same as the associated Kerr solution expansions writ- 
ten for the particle with parameters M, S. 

• If inequalities (|7j), (jHl) hold for the first particle and the parameters of 
the second particle satisfy inequalities 

^ ~ ^ ~ A^, (11) 

then in the PN approximation the second particle can be considered as a 
trial particle. In this case the presence of the second particle in the system 
impacts in no way on the dynamic equations for the first particle; as for 
the second particle, it moves in the field of gravity generated by the first 
particle like in the external field. In so doing the second particle is a trial 
spin particle. 

Input equations © are written in each order of approximation in the 
form of one equation 00, three equations Ok, and six equations mn. The 
consistency of arrangement of orders of smallness with equations (j21) 
shows up in the fact that written-out equation chain 



00; A 



2k 



Ok- A^'^+i 



mn: 



^2k+2 



allows us to sequentially determine 7 oo; ^ ofc, 7 mn- It is 



therewith assumed that coordinate conditions consistent with this procedure 
are used. 



5. Explicit form of the Kerr solution expansion 



Let us present the explicit form of the principal terms of the expansions of 
7q,/3 in smallness parameters u = m/r, v = a/r, where a = ^ (skSk)- The 
expansions are borrowed from ref. [20]. In the expressions for 700, 7ofc, 7mn 
we restrict ourselves to writing-out only those terms that contain multipliers 
given by u, v?. vu, vu^, v'^u, v^u. 

In the Kerr solution, the tensor of intrinsic angular momentum reduced 
to unit mass, Smn, has as few as one nonzero component S12, whose value 
coincides with a. Thus, 

S12 = a; S23 = S31 = 0; si = S2 = 0; S3 = a. (12) 

Axial momentum vector Sk is related to tensor Smn by the standard rela- 
tion 
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Smn - SmncSc, - ^E^abSab, S3 - U. (13) 

From (fT^ . (fT^ it follows that the following relations take place: 

1 = (14) 

cosH = ^'-^-l^^bXb) _ ^^^^ 
The expression for 700 in Cartesian coordinates is given by: 

,2 a SgSl, [ XgXb 1 



700 = 4m + - 6^ ( 2^ - ■v'-u + ... 
= 4^ + S - - \5ab) + ... . 



(16) 

The expressions for 701, 702, 7o3 in Cartesian coordinates are given by: 



7ofc = 2 



m{SkaXa) 



r-3 



r, m^{SkaXa) m{s kcXc) {SaXa){si,Xi,) , m{SkaXa) ( „ „ \ J_ 

(17) 

The expressions for 711, 722, 733,712, 713, 723 in Cartesian coordinates 
are given by: 



771 XfYiXji rt^^ r [(^Smc^c) •^n ~l~ (■^nc'^c) -^mj /io\ 
'ymn=^^ 2 — 5„„ + + .... (18) 



With vanishing Smn expressions flTH|) - (fTS|) transfer to Schwarzschild so- 
lution expansions. 

6. Construction of the Einstein equation solution cor- 
responding to two Kerr particles 

6.1. An approach to using the EIH method in the problem of two 
Kerr particles 

The approach differs basically in nothing from the one used for two parti- 
cles having nonzero mass, but possessing no intrinsic angular momentums 
(Schwarzschild particles). The problem specificity consists in the following: 
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First , in all orders of smallness functions 700, 
parts, 



7ofc, Imn split into two 



7oo 
7ofc 



7oo + 7oo 
7ofc + lok 



(19) 



The first parts (with circumflex accent) are solutions to the associated 
problem for two Schwarzschild particles. The second parts (with bar) are 
additions that owe their origin to the particle intrinsic momentums. Splitting 
(|19p can be performed always without loss of generality. 

Second , each of the parts of function (fT^ is expanded in smallness pa- 
rameters so that they coincide with the exact Kerr solution expansions in 
the two smallness parameters used in the EIH procedure. (Recall that in the 
problem of two Schwarzschild particles one smallness parameter was used.) 
The following arrangement of orders of smallness will be used according to 
exact solution expansions 



7oo = 
7ofc = 

7rran 



7oo + 7oo + 7oo + ••• 

.2 4 6 



+ 



7oo + 7oo + ••• 

L4 6 



7ofe + 7ofc + ••• 
.3 5 



+ 



7 ofc + 7 ofc + 
.3 5 



7mn + 7 
A 6 



+ 



7mn + 7 mn + 7mn + 
.4 5 6 



(20) 



In the expansion of for the exact solution the odd order appears in 
the approximation. 

Third , the exact solutions are written with the same coordinate conditions 
as in the EIH procedure. In our consideration these will be the harmonic 
coordinate conditions. 

Fourth , functions (pn|l are substituted both into the dynamic equations 
and coordinate conditions and the dynamic equations and coordinate condi- 
tions, which the second parts of the functions (with circumflex accent) have 
to satisfy, are determined. 

Fifth , the resultant equations and conditions are solved. In so doing the 
conditions of solvability of the dynamic equations and coordinate conditions 
are satisfied as well. 
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6.2. Determination of 7ojfc 

5 

The EIH procedures can be applied essentially in their standard form until 
7ofe has to be determined. We will not write out the associated expressions 

5 

as they are cumbersome. These expressions along with the details of their 
derivation will be presented in a preprint to be published. 

In this section we only present the explicit expression for 7ofc- The 

5 

function 7 ofc is determined from equation [Ok; X^] and coordinate condition 

5 

[c.c.;Qk; X^]. The calculation of surface integrals from the coordinate condi- 
tion leads to the following solvabihty conditions: 

M = -M {l;,^^ - -—- ,n = -,n(nn)--—. (21) 

As a result we arrive at: 

I K>0 ^1 -^i 

MHk I 2 "^^(^'-f^')^fc ^ mM{Xiii)Rk (22) 

R^n R^n R^n 

70. = 2^^^^ - 5 ^(SMis x.)is.x,) ^ ms^ ^g^g^^ 

5 K>0 1 1 ' 1 

M{S,X,)(Xii)Sk ^M(SA)Sk , ^M{Xiii){Skdi) 
+ t) -5 Z -g h Z -5 



rr r 



2, 



I A mM(skiRi) , mM(SkiRi) _ o ^l^'^O (^fe'^') ('231 

r, M^{SkiXi) _ mM{XiRi){SkiXi) , ^ mM{XaSaiRi,)Xu 
^ 7f ^Wrl ^W^^ 

I yt mMjskiXi) I -,r, mM{XaSabRb)Rk , a mM{XiRi){XaSabRb)Xk 

I f, mM{XaSabRb)Rk I /I mM{SkiXi) _ (. mM{XiRi){SkiRi) 
i?6ri "^^ RSri " RSri 

7. A method for derivation of the dynamic equations in 
the PN approximation 

The above expression for 7ofc is needed to derive the dynamic equations 

5 

determining the motion of centers of masses of particles as well as equations 
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determining time derivatives of spins, i.e. S^n a-nd |mn- 

The expression for 7ofc involves not S mn, s mm but S mm ^ mn- The equa- 

5 3333 

tions for S mn aiid s mn appear after 7 ofe has been substituted into the coor- 

33 5 

dinate condition relating 7 oit and 1 kssi and this condition is used to derive 

5 ' 6 ' 

1 mn- But we need not determine the complete expression for 1 mm it will 

6 6 

suffice to only determine the anti-symmetric part of 7 mn and require that it 

6 

vanish. It is this condition that will give the dynamic equations for the spin 
evolution. 

The equations of translational motion of particles are derived as before 
through integration of equation [mn; A^]. A new point is that, first, a different 
decomposition of the terms of this equation to the curl combination and the 

18 

part denoted byX) (a different amount of ttj and a different explicit form 

i=l 

18 

of them) takes place. Second, in the integration of ai it should be kept in 

j=i 

mind that the contributions are made not only by Schwarzschild parts '^fj,i,, 
but also by Kerr parts 7^j,. The strategy of derivation of the equations of 
translational motion of particles consists in determination of 7^,^, 7^,^, which 
make a contribution to the integral of [mn;A^], and in calculation of the 
integral. Note that the contribution of 7 mn to the integral is zero, so for this 

6 

purpose 7 mn need not be determined in the explicit form. 

6 

The resultant dynamic equations determining the motion of center of 
masses of particles coincide with those presented in many papers, e.g., [1]. 
Therefore, we do not write them out. Note only that the laws of conservation 
of energy and particle momentum vector follow from these equations. 



8. The dynamic equations for spins 

The resultant equations determining the spin dynamics that have been de- 
rived by the above method are given by: 

|c = 9^ [Riii) Sc + 2^ (Rim) - 2^ (Siii) R,. (24) 
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9. Time variation in total momentum 

The substitution of the rcsuhant dynamic equations for spins into the equa- 
tions for time derivative of total momentum of the particle system leads to 
the following expression for Mc'. 

Me = +^ {-2 (Siii) + 9 (Riii) + 2 (Rifji) S^} 

{-2 {sifji) R, + 9 {Rir]i) + 2 (i?,^) s^} 
+e.rr.n | _6 "^^^^-g^ ^ - A^^^^^R^ + 3^^i%^i?„} (26) 

f mM{S^iR,){R,Ri) mM{Smiii) „ ^ mM(S^;>?;) p ] 

-3^i{Si + si)Ri){Sa + s,i)Ri 

The analysis of the structure of the right-hand side of the equation for 
Mc suggests that no re- determination of the total angular momentum of the 
particle system can get the law of conservation of momentum to be obeyed 
in the PN approximation. 

For a spin particle system located in plane-asymptotic background space 
it may be possible to construct a conserved quantity similar to the total an- 
gular momentum. As it follows from our consideration, however, the quantity 
will not reduce definitely to the sum of orbital moments of the particles and 
their intrinsic angular momentums. 

The obtained result of nonconservation of the total angular momentum of 
the particle system in the PN approximation appears unusual, however, this 
result seems not to contradict the first principles of general relativity. Most 
probably, similar effects will appear with respect to energy and momentum 
as well at some approximation step. 

10. Numerical estimations for the gyroscopes on satel- 
lite GPB 

For the numerical estimations we will use the input data for gyroscopes that 

are presented in ref. [3]. According to that paper, the gyroscopes were 
spheres of radius Vg = 1.9cm and mass m = 75g. The initial angular velocity 
was a; = 27 000 rad/s. The satellite was launched onto a nearly ideal polar 
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orbit of radius R = 7027 km. The orbital velocity of the gyroscopes is 
V = 7.5 ■ 10^ cm/s. Ratio v/c is v/c = 2.5 ■ 10~^. The ratio of half the 
Schwarzschild radius of the Earth to the gyroscope orbit radius is 

_ = 6.5 ■ 10-^°. (27) 
R 

The gyroscope axes were in the satellite orbit plane. 

The equation used in ref. [3] to estimate the gyroscope axis of revolu- 
tion deviation is derived in that paper itself using Landau-Lifshitz energy- 
momentum pseudo-tensor [23]. It is given by: 

Sfc = {- (siRi) Rk - 19 (siRi) Rk + Wsk (RiRi) } . (28) 

Superpose the satellite orbit plane to plane {x, y) and write the change 
in radius-vector in form 

27r 

Ri = Rcos (nt) ; R2 = Rsin iyit) ; ^1 = ^- (29) 
From these conditions it follows that 

?7i = -RVtsin (nt) ; 7)2 = R^cos (nt) . (30) 

At the initial time the gyroscope axis of revolution orientation coincides 
with axis x, 

si{0) = So; 521(0) = 0. (31) 

From Anderson's equation it follows that the contribution to the time 
derivatives from the spin is made by the two first terms in the right-hand 
side of fj28p . Upon averaging of these terms over the orbit for the time of one 
revolution we obtain that 

Si S2 M 

—T = 0; — T = — . 32) 

So So R 

Thus, in one revolution around the Earth the gyroscope axes gets dis- 
placed by angle A9 of 

AO = 2M/R = 1.27 -lO^^. (33) 
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In its functioning the satellite GPB had to perform about 5000 revolu- 
tions, hence, the gyroscope deviation angles of the order of 6 ■ 10^^ radian 
had to be measured in the experiment. 

From equations (j^^ . (j^^ for spins that were derived with the method of 
this paper it follows that, first, the deviation must be 2 times as small and, 
second, the gyroscope axes deviate to the side opposite to the one following 
from equations (j2HI)- 

11. Conclusions 

As far as we know, the procedure of expansion of the Kerr solution written in 
the harmonic coordinates had been performed by nobody before our paper 
[20], hence, had not been used for normalization of expressions derived with 
the EIH method. This entitles us to suggest that the particle spin motion 
results presented elsewhere should be validated through comparison with the 
Kerr solution expansions in Section 5. 

Another test for correctness of the equations for spins that are presented 
in the literature is the absence of any anti-symmetric part in the expression 
for 7mn, which can appear from coordinate condition [c.c; A^]. 

6 

Now as for the amount of the gyroscope axis of revolution deviation on 
satellite GP-B. From the consideration presented in Section 10 it follows that 
our equations ((211), give /S.6 = —M/R for the deviation (per revolution 
around the Earth), while equations give A9 = +2M/R. We guess 
that the general relativity predicts the effect that agrees with our equations 
fl24|) . (j2Sl)- It is extremely interesting to learn the results of the forthcoming 
gyroscope experiment on satellite Gravity Probe B! 

From the obtained results it follows that the case in point should be not 
only the effects of the deviation from motion along the geodesic, precession 
and axis of revolution orientation variation, but also the change in lengths 
of intrinsic momentums, i.e. values of angular velocity of rotation. Gyro- 
scope rotation acceleration/deceleration on different segments of the world 
line could be noted in the laboratory without orbiting as well. To do this 
sufficient statistics should be compiled regarding the amount of revolutions 
made by the gyroscope about its axis for one and the same time on segments, 
where the spin projection onto the radius- vector is positive and negative. The 
comparison of the revolution amount on these two segments can serve a test 
for the general relativity validation. 
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